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Abstract 

The leptonic QED radiative corrections are calculated in the next-to-leading log approx¬ 
imation O [a 2 \n(Q 2 /mh)\ for unpolarized deeply inelastic ep-scattering in the case of 
mixed variables. The corrections are determined using mass factorization in the OMS- 
scherne for the double-differential scattering cross sections. 



1 Introduction 


Deeply inelastic electron-nucleon scattering provides one of the cleanest methods to test QCD 
by measuring the scaling violations of structure functions from which the QCD-scale Aqcd, 
one of the fundamental parameters of the Standard Model, is derived with high precision. The 
QCD-analysis also delivers the various twist-2 parton distributions as functions of x and Q 2 , 
the detailed knowledge of which is necessary for the forthcoming experimental search for the 
Higgs-boson and new, yet undiscovered states beyond the Standard Model, at TEVATRON 
and LHC. One of the major goals of the experiments HI and ZEUS at the ep collider HERA 
at DESY is to perform a QCD test at large space-like virtualities Q 2 at high precision. This 
presumes to know the QED radiative corrections to the double-differential scattering cross 
sections of deeply inelastic ep corrections as precisely as possible. Previous calculations of the 
radiative corrections for the unpolarized cross sections at leading order [l-5][], the leading-log 
level [7-12] to leading and higher orders, and QED-resummations of small-re terms [11,13]Q 
orders revealed that these corrections are very large, particularly in the interesting regions of 
small x and for almost all x at large inelasticities y = Q 2 /(Sx). The corrections are, moreover, 


complicated by a new type of sizeable contributions as the Compton-peak [|15 
necessary to extend the calculations to higher orders. 


This makes it 


The higher order leading-logarithmic contributions O ( aL) 


tained as the leading order solution of the associated renormalization group equations [16 


to QED corrections are ob- 

for 


mass factorization. These corrections are universal, process-independent w.r.t. their structure, 
and are given in terms of Mellin-convolutions of leading order QED splitting functions. The 
next-to-leading order (NLO) corrections can be obtained along the same line. However, be¬ 
sides the splitting functions to NLO also the respective process-dependent Wilson coefficients 
and operator matrix elements in the on-mass-shell (OMS) scheme contribute. In the past 
this method was applied to a single differential distribution for the 0(a 2 ) initial-state QED 
radiative corrections to e + e“ 


/j 1 1 i 


in Ref. [17 


In this letter we calculate the leptonic QED corrections to deeply inelastic ep scattering, 
where we define the double-differential scattering cross section choosing mixed variables [8,2,5], 
in 0(a 2 L). If compared to the case dealt with in [17[] to 0(a 2 L ) the present calculation is 
more complicated due to the emergence of final state radiation, the double-differential cross 
section and the relevant rescaling, which implies non-Mcllin type convolutions in general. After 
summarizing main kinematic aspects we present the different contributions to the leptonic 
NLO QED corrections. Further details of the calculation and numerical results will be given 
elsewhere (TB| . 


2 Mixed variables 

ep collider experiments allow to measure the kinematic variables defining the inclusive deep- 
inelastic scattering cross sections in various ways since four kinematic variables are available in 
principle with the energies and angles of both the outgoing lepton and the struck-quark, see 

1 For QED corrections to polarized lepton scattering off polarized nucleons see |J. 

2 For related numerical results for QCD singlet evolution see |fp|. 
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e.g. |T9|. At the Born level all methods are equivalent, however, resolution effects as a conse¬ 


quence of the detector’s structure, differ in certain kinematic regions. The Bremsstrahlung- 
effects of QED radiative corrections change this picture drastically and the QED correction 
factors depend on the way the kinematic variables, as e.g. Bjorken-?/ and the virtuality Q 2 
are measured.^ In the present paper we calculate the NLO radiative corrections in the case of 
neutral current deep-inelastic scattering for mixed variables, i.e. that Q 2 = Qf is measured at 
the leptonic and y = is measured at the hadronic vertex, and x m = Qf/(Syh)- The Born 
cross section for 7 -exchange is given by[] : 


with 


✓ A T (°) , N 2-ito 2 

da =C'(0>0)(j/,g 2 )= -.2 p t <a2 


dydQ 2 


yQ 3 4 


y 2xFi(x, Q ) + 2(1 - y) F 2 (x, Q ) 


( 1 ) 


Nf 
k =1 


Fl(x,Q 2 ) = [q k {x,Q 2 ) + q k (x,Q 2 


F 2 (x, Q 2 ) = 2xFi(x, Q 2 ) + F l (x, Q 2 ) . 


( 2 ) 

( 3 ) 


Here, Ei > 2 ,l(x, Q 2 ) denote the nucleon structure functions for photon exchange, and q(x,Q 2 ) 
and q(x, Q 2 ) are the quark- and antic^uark distribution functions. The sub-system variables 
obey the following rescaling relations for initial- and final-state radiation : 


ISR : y — Q 2 = zQ 2 = Ql, S = zS, x = zx m = x h , 

J\z) = 1, z£ = min |j/ h , ^|| , 

FSR : y = y h , Q 2 = = Q* § = S, x = ^ = x h , 

z z 

J F (z) = -, Zq = x m . 


( 4 ) 

( 5 ) 


Here, J I,F (z ) are the initial- and final-state Jacobians d 2 (y,Q 2 )/d 2 (yh,Q 2 ), and z 0 marks the 
lower bound of the sub-system rescaling variable z e [z 0 ,l]. The rescaling in Eqs. (^,|5|) was 
chosen such that both the initial- and final-state operator matrix elements can be expressed by 
a variable z e [0,1]. Ql is introduced as a scale to cut away contributions of the Compton peak. 
Although these terms do formally belong to the QED radiative corrections, they stem from a 
kinematic domain of low virtualities and are therefore not being associated to deep inelastic 
scattering. The scale Ql ~ 1 GeV 2 or larger can be chosen by experiment accepting only 
those events in the sample to be analyzed for which the hadronic Q 2 is larger than Ql. By this 
measure the Compton peak is cut away widely and the QED-correction factor is dominated 
by a deep-inelastic sub-process by far. Finally we note that in the case of mixed variables the 
leptonic QED radiative corrections can be easily grouped into those for the initial and final 
state. The separation scale between the two kinematic regions is Qf, i.e. we speak of initial 
state radiation for Ql < Ql and of final state radiation for Ql > Q 2 . 

3 See e.g. Refs. [8,2,5] for a comparison of a wide range of different choices of measurement. 

4 As we discuss radiative corrections we will speak about an 0(a k ) contribution for a cross section being 
proportional to a 2+k . 
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3 NLO corrections 


In this paper we limit the consideration to the calculation of the NLO corrections to leptonic 
variables for one-photon exchange in electron-nucleon scattering. This approach is widely 
model independent and allows to refer to general non-perturbative parameterizations of the 
structure functions which describe the hadronic tensor. In this way a direct unfolding of the 
experimentally measured structure functions is possible down to the range in Q 2 and x in 
which partonic approaches fail to provide a description of structure functions.^ The radiative 
corrections are thus valid as well for inclusive diffractive ep-scattering, see e.g. P2fl. The 


radiative corrections to be calculated concern a neutral-current process in which the lepton stays 
the same : both an incoming electron (positron) results into an outgoing electron (positron). 
As was already pointed out in Ref. |J changes of the fermion species are possible starting with 
2-loop order. To keep the QED A"-factor close to the basic process we limit the consideration of 
the present paper to those contributions in which the electron-type is preserved. The scattering 
cross section to 0(a 2 ) reads : 


d 2 


a 


dy h dQi 


(</-,«) + D {in ($) C^(y k ,Ql) + C^\y k ,Ql) 


( 6 ) 


+ (K ) 2 {hi 2 (§) (!/„,«?) + In (d) C^\y h ,Ql) + C^\y h ,Qt) 


with C l ' 0 ’ 0 \y h , Q 2 ) = d 2 a°/dyhdQ 2 and a = e 2 /(47r). The 0[(aL)\ and 0[(aL) 2 ] corrections 
were calculated in Ref. (T2|j . The term Qf) was derived in Ref. || completing the 

0(a) corrections. We re-calculated these corrections and agree with the previous results. 

The NLO-correction C^ 2,1 \y h , Q 2 ) can be obtained representing the scattering cross section 
using mass-factorization, see [|17]][]. Although the differential scattering cross section does not 
contain any mass singularity, one may decompose it in terms of Wilson coefficients and operator- 
matrix elements being convoluted with the Born cross section. In this decomposition both the 
operator matrix elements and the Wilson coefficients depend on the factorization scale /j 2 . One 
writes the scattering cross section as 


d 2 a 

dyhdQ 2 


d 2 a° 

dyhdQ 2 


r{ e ® d ee ® rf e + r 


I 

7e 



(7) 


with (z iff 2 j miff) the initial and final state operator matrix elements and <3y/(z, Q 2 / ff 2 ) the 

respective Wilson coefficients. ® denotes a convolution, which depends on specific rescalings 
of the chosen kinematic variables for the differential cross sections and is specified below. Both 
the operator matrix elements and the Wilson coefficients obey the representations 



( 8 ) 

( 9 ) 


5 Quarkonic and lepton-hadron interference contributions are not included in this approach. They are of 
smaller size in general |^o|, |j. ^l[ . 

6 Very recently also the electron energy spectrum in muon decay has been calculated using this method (23) . 
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The sequences {ij} and {kl} in the above do always denote j(l) for the incoming and i(k) the 
outgoing particle. As the differential cross section is /{-independent, the cross section is ex¬ 
pressed by convolutions of the functions T I i f'^ rn ' n \z) and a^ l,n \z) such that the /^-dependence 
cancels and a structure like Eq. (||) is obtained. The treatment in the OMS scheme assumes that 
the light fermion mass, m e , is kept everywhere it is giving a finite contribution to the scattering 
cross section if compared to the large scale Q 2 , i.e. the only terms being neglected are power 


corrections which are of O 


mrjQ 


2 \k 


k > 1 and therefore small. The last step is necessary 
to maintain the anticipated convolution structure which, in parts, is of the Mellin-type, as also 
in a massless approach. 

In the subsequent relations we make frequent use of the rescaling (fU|). For this purpose 
we introduce the following short-hand notation for the rescaling of a function F(y, Q 2 ) 


FjAv, Q 2 ) = F (y = y IiF , Q 2 = Q?, F ) , (10) 

where /, F label the respective type of rescaling. The NLO-corrections may be grouped into 
the following contributions : 


i LO initial and final state radiation off C ( F 1 \y^Q 2 ) 

ii coupling constant renormalization of C^l ,1 \y,Q 2 ) 

iii LO initial state splitting of P ie at C^\y, Q 2 ) 

iv LO final state splitting of P n at C^ 1 e ,1 \y,Q 2 ) 

v NLO initial and final state radiation off Cf e fi \y-> Q 2 ) 
The function C^ 2,1 \y ? Q 2 ) then is given by 


(ii) 

i=\ 

where y = yh, Q 2 = Q 2 . The function (y, Q 2 ) is obtained by 

c! v> (y,Q 2 ) = [ dzP°„(z) [0 {z - 4) J‘{z)C?» M 4 - C M(y, Q-) 

+ £ dzP e ° e (z) [6 (z - z£) J F (z)&4 1 \y,Q 2 )-C i - 1 0\y,Q 1 )], (12) 

where C^^^y, Q 2 ) denotes the non-logarithmic part of the 0(a) correction || |T8| and P f ° e (z) 
is the LO fermion-fermion splitting function 


po ( z ) = 1 + z " 

ee[) 1-Z 


To illustrate the convolution structure in (7) we note that 

cM( y ,Q 2 ) = |s<p> + +rfy-q ®c«»'°)(a.o 2 ). 


(13) 


(14) 
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Also the LO off-diagonal splitting functions 


p^z) = z 2 + (i-zy 

1 + (1 -zf 


PU Z ) = 


(15) 

(16) 


occur in other contributions to (7 2,1 ). Here both for LO and NLO splitting functions for equal 
particle transitions we write the contributions for z < 1 and account for the +-functions in 
explicit form below. 

We express the final result for C( 2,l \y, Q 2 ) in terms of the QED coupling constant a(m 2 ). 
Therefore we have to expand all contributions using 


a(fi 2 ) = a{m 2 ) 


l-ffthi 4 

47r V mi, 


( 17 ) 


with do = —4/3. In the present paper we only consider photonic corrections and those due to 
a single light fermion. The generalization to various fermions and non-perturbative contribu¬ 
tions due to light quarks is, however, straightforward. The running-coupling contribution to 
C {2A) (y,Q 2 ) resulting from C ( ' 11) (y, Q 2 ) is thus given by 


c 1 f' 1 ) fa.Q 2 ) = -§c' ( 1 ’ 1 | fa.e 2 ) 


(18) 


The contributions C^\y, Q 2 ) refer to two new 0(a) cross sections : d 2 a ie ^/dydQ 2 and 
^ 2 (7 e 7,(i) /dydQ 2 . The corrections are 


cf %,Q 2 ) = I dzP"(z)J‘(z)C l ‘-j’( V ,Q 2 ) , 


( 1 . 1 )/ 


7e\ 

•0 


e 7v 

dP’fe.Q 2 ) = [ r dzP%(z)j F (z)c { ;S(y,Q 2 ) 


(19) 

( 20 ) 


In these processes electrons and positrons are pair produced in the final state. The 0(a) 
sub-system cross sections read : 


dV 7e ’ (1) 

dydQ 2 

dV 7 ’ (1) 


a 
2 7T 
a 

27 


In (-yOf\y,Q*) + CV\y,Q->) 


In 


mt 


' Q 2 


mi 


Ciy o \y,Q 2 )+CilP(y,Q 


dydQ 2 

Here, the functions C-j ,0 \y,Q 2 ) are given by 

Ciy'fe.Q 2 ) = fdzP^(z)J,(z)cf fi \y,Q 2 ), 
Jz 0 

C% 0> (y,Q 2 ) = j'dzP° e (z)J F (z)c£'°\y,Q 2 ), 


( 21 ) 

( 22 ) 


(23) 

(24) 
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and the functions (y, Q 2 ) and CjF^(y, Q 2 ) are : 


C^\v,Q 2 ) = 


27TQ; 2 r 1 dz 

VhQi 1 4 z 2 


P°{z) In 


(! -z)(z~ y h ) 


(1 + 2z) In ( — 

\Vh, 


— 1 H-b 2^(1 — z ) 

Vh 

(z - yh)Pe 7 (z) In 

Vh 


+ 2 
+ (z- y h ) 


yfi ^XmFl^Xh, Qh) 


(! -z)(z- y h ) 


Vh In — 

,Vh, 


1 - 2y h + 82: + 2 y h z - 82:" 


F 2 (x h ,Q 2 h ) 


27ro; 2 f 1 


VhQ, 


dz\ 


1 Jz 0 


-<(*) In (V^) + f 1 + -) HVh) - - + — 

z 1 VI — z J \ z) z zy h 

X yh ‘^ , XmF\(Xh.i Qh) 


+2 


-P 7 e(z) (- - Vhj In (“y~7") + Vh In (y h ) + 


1 -Vh 2 1 ~ Vh 


c^M*) = 2 ^f ] d 4 


X F 2 (x h ,Ql)'j, 

1 


VhQ 


l Jz 0 


P"Jz) In 


z 

+ 2 
z - Vh 


1 -Vh 
1-z 


+ z In 



-1 + 


IJh 


(25) 


Vh C ^‘ X mF\{XhiQh) 


i^)(l-y fc )In(l-^)-y fc ln(A 

F 2 (x h ,Ql) 


27TQ: 2 f l , 

H-— / dz< 


VhQi Jzq 


2z 2 

Pfeiz) In 


(2 - y h - 3z + 2^2)) 


1, , N 0 3 1 

— In (yh) + 2 I 

2; 2; y h z 


X yh ‘^‘•^mFlip^hiQh) 
+ 2 


zPt( z )( l _ yh ) ln f —--— ] + 2/h ln (j/fc) 


:(1 - Vh){T~ yh ~ 142 + 4:y h z + 62 s - 2y h z 2 ) 


F 2 (x h ,Ql)\. 


(26) 


Finally, the term C^ 2,1 \y, Q 2 ) is given by : 


cf'W 2 ) = 


P 


l,JVS,OM 

ee,S 


(z) [0 (z - 4) J'(z)cr\y, Q 2 ) ~ C^°>(y, Q 2 ) 


p‘; p s s - OM (z)j'(z)cr\y,Q 2 ) 


R 


1, NS, OM 


ee,T 


M [» (2 - 2o F ) J F (z)C ( °-°\y,Q 2 ) - C<«>W) 




(0,0), 


(27) 
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Here P e } e $ T (z ) denote the space- and time-like NLO QED splitting functions of the non-singlet 
(NS) and pure-singlet (PS) channels. The label OM refers to the OMS scheme. We express 
these splitting functions referring to the splitting functions in the MS-scheme 


pi,ATS,MS 
z ee.S 


(z) = — 2P e ° e (z) ln(z) 


ln(l — z) + 


- \( l + z ) ln2 (^) - ^( 3 + 7 z) 111 ( 2 ) 


- 5(1 -*)-\\PLW 


ln (-) + 3j 


pi, PS, MS 
ee,S 


(z) = — (1 + 2 ) In 2 ( 2 ) + 
Plef m (z) = ^ e (z)\n(z) 


1 + 52 + - 2 2 


+ 2(1 — z 

2 / 1-2 


(28) 


ln(^) + - (—^) (10 + c + 28a 2 ) , (29) 


ln(l — z) — lll( 2 ) + - 


9 

+ i(l + 2 ) In 2 ( 2 ) - ^(7 + 32) 111 ( 2 ) 


-5(1 -z)--\P!e( z ) 


5 


111 ( 2 ) + - 


P e ™( 2 ) = (1 T 2 ) In 2 ( 2 ) 


5 T 92 T - 2 2 


+ 2 ( 1-2 
4/1-2 


111 ( 2 ) — 
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(30) 

) (5 + 23z + 14; 2 ) . (31) 


Note that the splitting functions P^sT^iz) refer to one light fermion or one light anti-fermion, 
contrary to the QCD singlet case, where always the sum of a light quark and an antiquark is 
considered. The splitting functions in the OMS are obtained by 


pi,TVS',OM/ 
*ee,S,T \ Z . 


pS,T,(U) 

ee 


M 


++/%-) + §W<w> w , 


1 + z * ( inf! 

1 ^ 

_ f - 2 V 

Z ^ 2/J 


(32) 

(33) 


with I = S,F = T, and P 


1, PS, OM 
ee,S,T 


» = nS MS ( 


4 Conclusions 

We calculated the 0(a 2 L) leptonic QED corrections to deep inelastic ep scattering for the 
case of mixed variables. The corrections are given in terms of double-differential distributions 
to be compared to the double differential Born cross section. The calculation was performed 
using the renormalization-group decomposition of the 2 -loop corrections to the differential 
cross section w.r.t. mass factorization in the OMS scheme for the light fermion mass. Mass 
factorization introduces an artificial factorization scale p 2 on which the physical cross section 
does not depend. By the elimination of this scale the structure of the differential cross section 
is re-organized in terms of pieces which can be calculated first individually and which are then 
composed to the desired correction term to be obtained. The NLO correction consists out of 
five terms, the LO ISR and FSR radiation correction of the non-logarithmic 0(a) contribution, 
a respective term due to charge renormalization, two new terms containing e — 7 — e initial and 
final state transitions, and the ISR and FSR OMS-NLO radiation correction to the Born term. 
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